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Abstract 

We study the electromagnetic field equations on an arbitrary quan- 
tum curved background in the semiclassical approximation of Loop 
Quantum Gravity. The effective interaction hamiltonian for the Maxwell 
and gravitational fields is obtained and the corresponding field equa- 
tions, which can be expressed as a modified wave equation for the 
Maxwell potential, are derived. We use these results to analyze elec- 
tromagnetic wave propagation on a quantum Robertson- Walker space 
time and show that Lorentz Invariance is not preserved. The formal- 
ism developed can be applied to the case where back reaction effects 
on the metric due to the electromagnetic field are taken into account, 
leading to non covariant field equations. 



1 Introduction 



There are several claims in the literature that one of the milestones of physics 
in the past century, namely Lorentz invariance, will no longer be true once a 
reliable quantum theory of gravity is achieved. 

In its most simple and geometrical form, this invariance can be realized 
by introducing the Minkowski metric rj^ as the "fixed backstage" to define 
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distances or norms in special relativity. For example a photon is a particle 
whose 4-velocity k^ has zero norm, i.e. 

Vfiuk^k = 0, 

and the above equation is the geometrical casting of the invariance of the 
speed of light. What if light travels through cosmological distances? No 
problem, according to GR we simply replace the flat metric for the nontrivial 
metric g^ v of the curved spacetime. Thus, Lorentz invariance in a local frame 
is generalized to covariance of the equations in an arbitrary frame. From this 
point of view the above equation is a scalar with respect to any coordinate 
transformation. 

What if the classical metric g^ is replaced by a quantum operator g^ u 7 
Here we face a problem, namely, the meaning of the resulting equation. A 
scalar equation can be obtained by taking an expectation value of the oper- 
ator g^ u with a suitable gravity state. If k^ is a null vector with respect to 
(g^u) for a particular state then in general it will not be a null vector with 
respect to a different gravity state. However, in this case one is not really 
breaking Lorentz invariance since taking expectation values with different 
states is like having classical metrics with different conformal structure and 
a vector is null with respect to both metrics if and only if the conformal 
structures are the same. 

As a matter of fact it is difficult to imagine how to obtain non covariant 
observables if one is able to construct a theory with covariant field equa- 
tions for the operators and the state vectors are invariant under the gauge 
transformations resulting from the diffeomorphism group. 

However, to this day we do not have such a theory. Rather, the leading 
candidates for a quantum theory of gravity provide models for light propaga- 
tion on a geometry constructed from a semiclassical quantum gravity state 
that break Lorentz invariance [U [2]. The effects predicted by these mod- 
els are within the detection limits of present technologies but the observed 
data severely compromise the validity of their results. Both loop quantum 
gravity and superstring theory predict a frequency dependent speed for pho- 
tons propagating on a quantum spacetime. However, the approximately 3000 
GRBs observed by BATSE and other space instruments showed that all the 
photons emitted by these bursts have the same flight time within the in- 
struments sensitivity limits. Another prediction of loop quantum gravity 
concerns the rotation of the polarization direction for linearly polarized ra- 
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diation but the observed synchrotron radiation for sources located at cos- 
mological distances put severe restrictions on the phenomenological coupling 
constant of the model [3l H] . It appears very likely that Lorentz invariance is 
preserved at the linearized approximation since the observed evidence points 
in this direction. 

Upon a close look at the two models one finds possible resolutions to the 
conundrums. For example, the dispersion relation for photons obtained from 
the superstring model can be written as [5] 

g^(k)W = 0, 

that is, the metric of the spacetime depends on the energy of the photon. 
The equation is fully covariant and it does not say that different photons 
move with different speeds. It only says that, for the specific spacetime con- 
structed with one photon interacting with gravity, the metric also depends 
on the energy of the photon. If one changes the energy of the photon one is 
also changing the spacetime and thus, it is "illegal" in GR to compare results 
coming from two different spacetimes. One could try to solve the problem of 
two photons with different energies interacting with a quantum gravity state 
but it appears very likely that if the problem is well set both photons should 
follow null geodesies as the observations suggest. 

The loop quantum gravity model also admits a second look that offers 
an alternative explanation for the unobserved prediction. If we assume that 
the classical electromagnetic field is the 2-form F^, Lorentz invariance is 
preserved [6], whereas if we assume the canonical variables (A^, —E u ) can be 
regarded as classical objects, the invariance is broken. The problem arises 
from the relationship 

When the metric becomes a quantum operator 'g fiu and one takes an expec- 
tation value of the above equation with a pure gravity state then both the 
covariant and contravariant versions of the electric field cannot be regarded 
as classical variables. If one assumes that the field of the l.h.s. is not affected 
by the gravity state one obtains non-lorentzian equations of motion. If, one 
the other hand, one assumes that the electric field of the r.h.s is transparent 
to the action of the expectation value then Lorentz invariance is preserved 
at a semiclassical approximation, i.e., neglecting back reaction effects. If the 
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invariance is broken it must come from taking into account these back reac- 
tion effects. 

The aim of this work is to study the propagation of light interacting 
with a semiclassical quantum gravity state when back reaction effects are 
included in this interaction. The goal is to see whether or not covariant 
equations of motion are obtained for this propagation. A more precise or 
technical meaning of this problem is given in Section 3 were it is defined and 
analyzed. However, in section 2 we address the propagation of photons on 
a non-flat semiclassical gravity state since these results are later used in the 
main section of this work. The derivations on Section 2 can also be used as 
a review of our previous work or as a toy model for photons propagating on 
a geometry given by semiclassical quantum gravity states that are peaked 
around a non trivial classical metric, as for example a Robertson Walker 
spacetime. Finally, in the Conclusions we summarize our work. 

2 Propagation of electromagnetic radiation 
on a non-flat semiclassical geometry. 

In this section we analyze the interaction of quantum gravity and Maxwell 
fields acting on quantum states that are a direct product of coherent states 
for the electromagnetic field and weave states for gravity. The non-trivial 
equations of motion that are obtained from such a scheme is called the semi- 
classical approximation of loop quantum gravity. A derivation of these equa- 
tions follows 

2.1 The effective interaction Hamiltonian in the semi- 
classical approximation 

Assume an arbitrary background metric g^ u , and consider a 3+1 splitting of 
the spacetime by introducing a foliation of space-like hypersurfaces. We can 
set coordinates (t, x) adapted to the foliation such that the lapse and shift 
N and iV a are 1 and respectively (particular gauge choice). Let q a b be the 
induced 3-metric on the S t (corresponding to t = const) surface^. 

1 From now on we will use Latin indices to denote spacial components and Greek indices 
for 4-dimensional components. 
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We will assume that there exists a geometric weave state |A) on S t such 
that, given the classical metric q a t>, (A\q ab \A) = q ab + where q ab is the 
quantum operator associated to the metric tensor, £p is Planck's length, and 
A is the typical length of the weave |A). Such a state could be constructed, 
for example, by introducing random oriented Planck scale circular loops that 
form a graph adapted to the local geometry, and considering the product of 
the traces of the holonomies along these loops [7]. 

Now, if E a and B a are the electric and magnetic (purely spacial) fields on 
that background, the Hamiltonian density that couples these fields to gravity 
is given by 

H EB = -/ d?xq ah ^det(q) (E a E b + B a B b ) 

= \ \ d 3 xq^(e a e b + b a b b ), (1) 

where in the last line we have rewritten the Hamiltonian in terms of the 
vector densities e a and b a associated with E a and B a , and q^ is the 3-metric 

divided by its determinant. When the Hamiltonian is expressed in these 
variables, it is possible to implement Thiemann's regularization procedure, 
which consists of a point splitting method where the operator associated 
with the metric divided by its determinant is written as the product of two 
operators w^(x) (each one given by the commutator of the Ashtekar connec- 
tion A l a and the square root of the volume operator associated with q ab , i.e. 

Wa(x) = A l a (x), \/V{x) [51 [9]) evaluated at different points. In this way, 
the quantum operator corresponding to the electric part of the Hamiltonian 
density is 

H E = \J d*xj Sy5 u wl{x)w[{y)e\x)e b {^Ux - y), (2) 

and similar for the magnetic part. Here, f t (x — y) is a regularization 
function that tends to 5(x — y) as e — > 0. The next step in the regulariza- 
tion procedure consists of introducing a triangulation of the hypersurface S t 
into tetrahedra adapted to the graph associated with the weave state |A) 
considered [SI E] • 
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The effective interaction Hamiltonian is then defined as the expectation 
value of the above operator in a semiclassical state that is given by the weave 
described before for the gravitational sector, and we will assume, in addition, 
that this state is close to a coherent state for the Maxwell sector, in such a way 
that, within our approximation, we can consider the electromagnetic field as 
a classical quantity. Under these assumptions, the effective Hamiltonian is 
given by 

He = iy"d 3 xy"A^/(A|u;^f)u;i(y)|A)e a (f)e b (y)/ e (f-y) 

= ^^(A|^(^)^(^)|A)e a (^)e 6 K), (3) 

Vi,Vj 

since, by construction, the operators w*(x) only act at the vertices Vi of 
the graph. 

If we now assume that the variation scale of e a is large compared to the 
typical length of the weave state, we can expand it in a Taylor series around 
the central point x of the graph. Keeping only linear terms in £p, and using 
the fact that the Hamiltonian is invariant under spacial rotations, we find 
that the electric part of the effective Hamiltonian is given, up to linear order, 
by: 

H E eff = Vdet(q) (^q ab E a E b + ^ py /det(q)e abc E a V b E^ , (4) 

and similar for B a . Here £ is a phenomenological coupling constant, e abc 
is the total antisymmetric tensor that represents the volume element of the 
classic 3-metric q a b, V& is the 3-dimensional covariant derivative consistent 
with q a b and we have expressed the Hamiltonian in the original variables 
(instead of vector densities). Also, we used the form E a and not the vector 
E a for reasons that will become clear in the following section. 

2.2 The field equations 

To derive the field equations from the above Hamiltonian we first need to 
determine the relationship between the electric and magnetic fields and the 
canonical variables. Let A be the Maxwell connection and F = dA the 
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associated Maxwell 2-form. Then, the electric and magnetic fields that an 
observer with 4-velocity t^ 1 would measure are given by 



Efx = F^t\ (5) 
B, = - X -e,fF p& t\ (6) 

respectively, where e^ up s is the totally antisymmetric tensor associated to 
the volume element of g^ v . 

Note that, form the above expressions, both E p and are (as stated be- 
fore) purely spacial vectors and that the definition ([5]) says that E a = —d t A a 
(we have chosen A = since we are only interested in wave propagation), 
which is independent of any background metric. On the other hand, if n a is 
the canonical momentum conjugated to A a , Hamilton equations are 

piTT 

- ^ (7) 
These expressions are consistent with (jSJ) only if 

dH = dHm 

which can be seen as a differential equation for E a (n b ). By solving it 
we obtain the following relation between the electric field and the canonical 
momentum 

ir a = -y/ddfijj (q ah E b + 2^£p^t(qj e abc V c £ b ) , (10) 
or, symbolically 

n a = -H ab E b , (11) 
where the "metric" operator H ab is given by 

H ab = ^Jdet{q) (q ab + 2^ P y/det(q) e abc V c ) . (12) 



Note that eq. ( fill) is the semiclassical analogue of the relationship ir a 
■^det(q)q ab E b . 
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Now, the second Hamilton equation, ([7]), leads to 

d t (H ab E b ) = H ab e b cd V c B d . (13) 

Introducing the definitions (jSJ) and ([6]) of E a and B a in terms of the 
Maxwell potential we obtain the following expression 

WF^ = 2^ P d t [det(q) 2 ]H ad e bcd V b (F 0c ), (14) 

where V M is the 4-dimensional covariant derivative associated with the 
background metric g and H ab is the inverse of H ab . From this equation it 
becomes clear that the term on the left is the spatial component of a covariant 
expression and the term that breaks covariance appears in the right hand side. 
In particular, notice that for a stationary metric this term vanishes and we 
obtain a Lorentz invariant propagation. 

2.3 Light propagation on a quantum flat FRW back- 
ground 

We can apply, as an example, the formalism developed in the previous section 
to the particular case of a flat Friedman- Robert son- Walker spacetime: 

ds 2 = -dt 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ). (15) 
In this case, the metric operator acting on a vector C a is given by 

H ab C b = (a{t)5 ab + 2^£ P a{t)h abc d c ) C b (16) 

with e abc the Levi-Civita symbol, and the corresponding field equation 
([IB"]) in the semiclassical approximation is 

dt(a(t)i) = VxS, (17) 

where we have adopted, for simplicity, vectorial notation and, for any 
given vector C, the quantity C is defined as: 

6 = C - 2^ P a{t) 2 V x C, (18) 

with V x C the usual 3-D curl in flat coordinates, i.e. x C\ = e bc d b C c . 
This is the only non trivial equation, since the other Hamilton equation (??) 
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gives no new information, it is just the definition of the electric field in terms 
of the potential. Note that in the above expression the indices are raised 
and lowered with the 5 a b 3-metric, and all the time dependence has been put 
explicitly. 

On the other hand, the classical counterpart of eq (TTTT) is 

d t (a(t)E) = VxB. (19) 

There is another equation that relates E and B, and it comes form the 
fact that, since F = dA, dF = 0, which written in terms of the electric and 
magnetic fields defined by ([5]) and ([6]) reads: 

a - V x E 

d t B = — B 20 

a a 

where a = This expression holds both classically and in the semi clas- 
sical approximation we are analyzing. 

Moreover, from the field equation (fTTl) we can derive the wave equation 
satisfied by the Maxwell potential: 

d t (ad t A) - ^ = 2&p (<9 t (a 3 <9 t (V x A)) - aV x (V 2 i)) , (21) 
while the classical covariant version is given by 

d t (ad t A) - — = 0. (22) 



2.3.1 Plane Wave solutions 

Consider, as a first step, light coming from a sufficiently close source. In that 
case we can take a(t) ~ const and the wave equation (J2TJ) reduces to 

a$A -YlA = 2££ P a 2 V x [ ad? A - — | . (23) 

If we propose a plane wave solution of the form A = Re (^A e lS J and in- 
troduce it in the above expression, we obtain the following dispersion relation 
for the wave vector k ll = 5 M 
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{-a 2 k 2 + k 2 ){l-2^ P a 2 k) = 0, (24) 

where k 2 = k 2 + k 2 + k 2 . We see, then, that the wave vector k^ satisfies the 
usual dispersion relation g^k^kv = 0. 

Note that this result holds in the geometric approximation of wave prop- 
agation, since in that case we are dealing with the high frequency limit and 
we can neglect time derivatives of the expansion factor compared to time 
derivatives of the electromagnetic field. 

If, on the other hand, we consider a source located at a cosmological dis- 
tance, then we can not take a(t) ~ const, we must take into account the 
terms containing time derivatives of a and solve the complete wave equa- 
tion f[2"Tj) perturbatively since it does nos admit plane wave solutions with 
constant amplitude of the form proposed above. This is done more easily if 
we introduce the conformal time rj, such that -5 = a -1 . Expressed in this 
conformal time, the FRW line element reads 

ds 2 = a(r]) 2 (-dv 2 + dx 2 + dy 2 + dz 2 ), (25) 
and the wave equation (!2~T!) can be rewritten as 

(l-2tl£pa 2 Vx)nA = AtlipaaV x A', (26) 



where OA = d 2 A — V 2 A and prime denotes derivative with respect to the 
conformal time rj. Note that, in view of the conformal flatness of the metric 
(|25|) . the associated classical wave equation is simply OA = 0. 

We will try to solve eq. (126]) in a perturbative way, by proposing a solution 
of the form 

A = A dass + C£pA, (27) 

where A dass is the classical plane wave solution A dass = Re(A e^ UJV ~ k x ^), 
with oj 2 = k 2 . Introducing the solution fTSTl) into (126]) and dropping terms of 

order (£^p) 2 we obtain the equation for A: 

OA = 4aa'V x A' dass 

= -Ataa'k x A' dass (28) 
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To find the final solution, we must say something about the expansion factor 
a(t). Consider an Einstein-De Sitter model, that is, FRW universe dominated 
by matter. In this case the radius of the universe is given by a(t) = at 2 ^ 3 , 
which, written in terms of the conformal time r\ is 

a(r ] ) = a(^(r ] -r ) o) + tl /3 y , (29) 

where the subindex denotes the moment to of emission of the elec- 
tromagnetic radiation, at which the expansion factor is assumed to have the 
value ao (this means that a = aot 2 ^ 3 ). Inserting this into (|28|) we can obtain 
the final solution: 



A = Re 
= Re 



(1 + ^p« 2 (» x A )(L( W , r) - L( W , r )))] , (30) 
l ^- l£ \A + C£pA(v^))} , (31) 



where ft = t|t and the function L(u, r) is given by 

a 3 r 2ar 3 . ( a 2 r 2 



^•''» = -^ + ^ +! (i^-" r4 J' (32) 

with r = f (?7 — ?7o) + ^c/ 3 o r ; i n terms of the comoving time, r = t 1 / 3 . 
Note that we have imposed the condition that the solution coincides with 
the classical wave at the emission time to. 

We see from eq. (I3"U1) that the final solution corresponds to a plane wave 
with the standard dispersion relation and hence with no modification on 
the propagation speed, but with a corrected amplitude vector due to quan- 
tum gravity effects. As we will see in the following, the interaction of the 
electromagnetic radiation with the quantum spacetime induces a frequency 
dependent correction on the polarization direction of the initial wave, while 
its amplitude is, within the linear approximation considered, not affected. 

To see how this correction behaves, consider a normalized and linearly 
polarized initial wave. Introducing a unitary right handed basis (ei,e2,n) 
and assuming that the initial Maxwell potential is polarized along the e\ 
direction, we obtain, by taking the real part of f[3"0"j) 
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A 



[e-i + t;£pe2Re(L(uj, r) — L(u, r ))]cos(ur] — k ■ x) 
— £,£peia 2 Im(L(u, r) — L(u, r ))sin(ur] — k ■ x). 



(33) 



From this expression we can not conclude any quantitative results for the 
correction to the magnitude of the amplitude vector, since it is given by 



and we have been dealing with the linear approximation only (and hence 
dropping quadratic terms in the whole calculation that lead to this equation). 
On the other hand, the angle of rotation 9 of the polarization vector can be 
obtained, leading to 



tan(0) = ZipaotQ 2/3 y/lm(L(u, r) - L(u, r )) 2 + Re(L(u, r) - L(u, r )) 2 . 



Here 6 is measured from the initial polarization direction, and it is an in- 
creasing function of both t and u. More precisely, from (|35[) it is possible to 
prove that the tangent of the polarization angle has, for a given frequency, 
a behavior where the dominant term is of the form t 4//3 (see fig [TJ), while for 
a fixed instant of time the angle grows linearly with the photon energy (fig 
12]). Notice that this result is different from the one obtained by Gambini and 
Pullin [2], where the dependence of the polarization angle was quadratic in 
the photon energy. Note, in addition, that the solution ( |30l) does not show 
the birefringence effect predicted in [2], since in our formalism the propaga- 
tion velocity does not depend on the frequency, nor on the polarization state 
of the wave. 



\A + Z£ P A(r},u)\ = \A \ + Om P ) 2 ) 



(34) 



(35) 
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Figure 1: Tangent of the polarization angle as a function of the comoving time 
for a given photon energy, the behavior is of the form tan(0) oc t 4 / 3 . 
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Figure 2: Tangent of the polarization angle as a function of the wave frequency 
for fixed t. The graph shows a linear behavior. 
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From the above considerations, we conclude that, through its flight time, 
the polarization direction of an electromagnetic wave will be rotating a fre- 
quency dependent angle 9{u). If a source emits a wave packet with a con- 
tinuum spectrum, the high frequency photons will rotate a larger angle than 
the less energetic ones, the net result being a loss of linear polarization. The 
fact that we observe, nevertheless, light coming from cosmological sources 
with a high level of linear polarization is indicative that the effect, if present, 
is very small (some orders of magnitude below the sensitivity of the current 
instruments). Even more, we could use the available observational data to 
put a bound on the coupling constant £. Clearly, that bound would differ 
form the value obtained in [21 H] , since there a quadratic effect was assumed. 

Just for completeness, suppose two photons are emitted simultaneously 
by a cosmological source located at distance L with identical (linear) polar- 
ization state, and with wavelengths Ai and A2. Then, at the detection time 
their respective polarization directions would rotate in such a way that the 
difference between the corresponding angles would be given by 



To derive this expression we used (1331) and, based on observational argu- 
ments, assumed ^ < 1 and therefore tan(#) ~ 9. We have also considered 
only the leading terms in (I35I) since we are dealing with high frequencies and 
large distances. Hence, by using (136|) we could, if we knew the constant a, 
obtain the desired bound for £. However, we do not have enough information 
on the parameters of the universe to determine a, which makes it extremely 
difficult to estimate that bound. On the other hand, for the curvature ef- 
fects to be appreciable, we should have observations of sources at very high 
redshift, and we believe that, for the available data, the flat background ap- 
proximation suffices and the bound obtained in [3], H] is the most reliable. 

Remarks: 

• One thing worth mentioning is that, for non-stationary geometries, 
the assumption that it is possible to construct a weave state that is 
peaked at that specific metric for all times is a very strong one. If 




(36) 
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the 3-metric has, indeed, a non-trivial evolution, it is encoded in its 
conjugated variable (directly related to the extrinsic curvature of the 
hypersurface) and, since the weave does not satisfy the properties of 
a coherent sate (namely, that it approximates the configuration vari- 
able and its canonical conjugate), when one takes expectation value 
there is no control on d t (q a b) and it could happen that (g a &) deviates 
significantly form the classical value after a short time. 

• Even more, there are indications that for arbitrary curved spaces, the 
weave states might not be solutions of the Hamiltonian Constraint 
[7]. Thus, if the weave states cannot be considered physical, the re- 
sults showed in this section may not reflect the real propagation of 
light on a semiclassical FRW spacetime, assuming of course one can 
find a suitable definition of semiclassical states. 

• Note also that, in the derivation of the above equations, we have not 
taken into account the back reaction effects on the background ge- 
ometry. It could happen that, even if the resulting equations locally 
preserve Lorentz invariance (for example if one only considers station- 
ary backgrounds), this invariance could be broken if one considers the 
back reaction effect. This possibility is analyzed in the next section. 

3 Lorentz Invariance and back reaction ef- 
fects 

In the previous sections we have considered wave propagation on a fixed 
background geometry, that is, we have neglected back reaction effects on the 
metric due to the electromagnetic field. However, Einstein's equations couple 
gravity to any other forms of energy, in particular, with the Maxwell field, 
and we expect, therefore, that not only the quantum geometry will affect the 
propagation of electromagnetic waves, but also that the latter will modify, in 
turn, the spacetime itself. Here we will try to account for this back reaction 
effect by applying the formalism developed in section 3 to the full Einstein- 
Maxwell theory. To do so, we first have to analyze if it is reasonable to 
consider that we are within the assumptions made in that section, namely, 
that we have a classical metric expressed in the appropriate gauge, and a 
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weave state that approximates that particular metric. 



The classical equations that describe the Einstein-Maxwell theory are 

Gfiv = SwT^, (37) 
Vfjm = 0, (38) 

where the Einstein tensor is determined by the stess-energy tensor 
associated to the Maxwell field F^, that is 

T, v = ^ (f»„F u ° - ±g^F 5p F 5 ^ , (39) 

and V M is the covariant derivative consistent with g pv . 

The idea is to solve this system of equations with an approximation 
scheme were at the zeroth order the electromagnetic field propagates on a 
flat background. In other words the deviation from flatness arises from the 
electromagnetic stress energy tensor. 

As stated above, there are two things that need to be considered: 

i) that there is a classical metric solution of f[3"Tj) written in the gauge 
where g tt = —1 and g ta = (that is, we introduce a foliation adapted 
to the time field £ M , such that the lapse and shift functions are N = 1 
and iV a = respectively), and 

ii) that it is reasonable to assume that we can construct a semiclassi- 
cal state \vf)) that approximates the 3-metric q a b induced in the spa- 
cial hypersurfaces of the above mentioned foliation, that is, such that 
{ip\q a b\^) = (lab plus corrections of order tp. 

If we are able to affirm these two statements, then we can assume that the 
field equation for the electromagnetic field, within the semiclassical approxi- 
mation of LQG, is given by eq. (TT3l) . with H ab and e a bc the ones corresponding 
to the 3-metric q a b associated to the solution of ( 1371) . 

The classical metric 

We have to solve the field equation fl37j) . with T^ y given by fl39l) . Since 
by assumption the Maxwell field introduces a small perturbation to the 
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Minkowski metric rj^, we will assume that the electromagnetic tensor is 
given by eF^, with e a small parameter that will allow us to solve the equa- 
tions by means of a perturbation expansion. Then, the classical metric will 
be given by 

g^u = v + + e2 a$ + ■■■ (40) 

Besides, we will set the data (in the initial surface) of gfi) , gj?J, etc, equal 
to zero, since we are interested in the case where there are no incoming 
gravitational waves, and all the perturbation is generated by the interaction 
with the Maxwell field. 

If we introduce (1401) into (1370 . and consider a similar perturbative ex- 
pansion for the electromagnetic field (and hence for the stress-energy tensor) 
given by 



F^ = e(F«+eF« + ...), (41) 
T„ v = e\T$+eT$ + ...), (42) 

we can solve the coupled equations order by order in a recursive way, 
such that each order is determined by the previous ones. In this work we will 
focus on the first non trivial corrections to the free fields. 

The term of Equation (13 7p that corresponds to linear order in e is simply 



Gg> = SttTW (43) 
= 0, (44) 

since, the first non zero term of the stress-energy tensor is of order e 2 , 
and, hence, we get the familiar result form linearized gravity [UJ 

d 8 d { ,g% - \ d 5 d s g$ - \ v d<>d s gV = 0, (45) 

where g$ = g$ — \v^g^ with = rj^g^J. Of course the above 
equation reduces to the well known wave equation cftdsgffl = in the Lorentz 
gauge. We will not, however, consider that gauge but another one consistent 
with ( fl3l) . We will go back to this issue later. 



17 



The only solution of eq. gSD is g$ = 0, since, as mentioned before, the 
data for this perturbation in the initial surface is zero. On the other hand, 
looking at the first order in e of eq (1551) we get 

= 0, (46) 

which sates, as expected, that Ffi) is the free Maxwell field in flat back- 
ground. 

Hence, the non trivial correction for both the electromagnetic field and 
metric tensor are at least of order e 2 . The field equation for gffi is, in analogy 

to (USD, 

&d { »9% ~ \ d 5 dsg$ - \ dotfgW = 8nT$ , (47) 

with 

T<8 = I (f$FU° - \v^F^ , (48) 

and where all the indices are raised and lowered with rj^. The above 
expressions prove that the first correction to the metric due to back reaction 
effects is generated by the free Maxwell filed. We will write, for simplicity, 

5V = V + e2 7/^ (49) 

with 7^ such that 7^ is a solution of ([47] ) . 

Now, the formalism developed in the previous section requires that the 
metric is written in coordinates in which it takes the form 

ds 2 = -dt 2 + q ab dx a dx b . (50) 

This can be easily done if we use the gauge freedom in eq (T47]) that 
corresponds, precisely, to a coordinate choice. It is well known that, under 
the transformation a: M 1— > x M + e the perturbation changes according to 
Jiw l— > Jfiu + dfi£v + dv£n- in order to satisfy (150]) . the gauge choice must 
be such that j tu = 0, which gives the four necessary conditions to fix the 
transformation generator namely 
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dttv + Quit = -Itu- (51) 

Hence, it is always possible to express the metric tensor in the form (|50p , 
and, in that particular gauge, 7^ has only spacial components. 

To obtain the corrections to the electromagnetic field due to back reaction 
effects, we insert the solution (j4"9~|) in the Maxwell equation V^F^ V = 0. By 
doing so, we obtain 

PFjS = 0, (52) 
W (3) = Y"d.FU+rr(T$'FV+r$'FV), (53) 

with the Christoffel symbol given by 

rg* = + d ul , 5 - d 5lllv ). (54) 

The above equations tell us that the first non-trivial correction to the 

(3) 

Maxwell field is F^ u , since we can use an analogous argument to the one 
used to conclude that g$ = 0; namely, we consider that the only incoming 
waves are the ones given by the free Maxwell field F$ and hence, since there 
is no source to generate it, F^J must be zero. Therefore, the electric and 
magnetic fields will be of the form 



E a = eE^ + e 3 E^ + ..., (55) 
B a = eB^ + e"B^ + .... (56) 

In the perturbative formalism described, we saw that each order in the 
perturbative expansion can be obtained from the previous ones, both for the 
metric and for the electromagnetic field. We will stop the calculation here, 
however, because we are just interested in the first non-trivial corrections, 
which, for the Einstein-Maxwell theory, are given by eqs fj47|) and fl52|) . 

We have proved statement i), namely, that it is possible to find a classi- 
cal metric solution of Einstein-Maxwell equations and write it in the gauge 
where g a = —1 and g ta = 0. It only remains to prove the existence of a 
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semiclassical state that approximates that metric (statement ii). 



The semiclassical state 

We assume there exists a semiclassical state that satisfies the following 
condition: 

• It is peaked at the classical 3- metric, i.e., 

{lp\%b\lp) = qab class + 0(£ P ) (57) 

= 6 ab + e 2 7ab + 0(£ P ). (58) 
We will derive here a formal solution of the above stated condition. 
In the perturbative approach we are considering, the Hamiltonian is given 

by 

H = H + e 2 H 2 , (59) 

where T~Lq is the Hamiltonian constraint corresponding to pure gravity, 
and TC 2 is the perturbation introduced by the Maxwell field (which is the 
usual u q a b(E a E b + B a B b ) n term) and that describes the coupling between 
the two fields. Similarly, we propose a semiclassical state of the form 

|^) = |^ )+e 2 |^ 2 ), (60) 

where j^o) is the known state for the unperturbed Hamiltonian (that 
is, for instance, the weave state |A)), and hence satisfies Ho\^o) = and 
{^o\q a b\^o) = S a b + 0(£p), while \ip 2 ) is a correction due to the electromag- 
netic field. 

Now, in order that the peakedness condition stated above be satisfied, 
the perturbation \ip 2 ) must be such that 

2Re{(^\q ab \^ 2 )) = lab (61) 

modulo corrections of the order of £p. 

In the following we will assume that it is possible to find a state \ip 2 ) 
that satisfies this expression and, hence, approximates the classical metric 
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derived in the previous subsection, up to corrections of order £p. Hence, we 
are within the assumptions made to derive eq. ([TBI) , which allows us to apply 
the formalism developed in section 3. 

3.1 The field equations: semiclassical photon propaga- 
tion with back reaction 

Recall from section 3 that the field equations are given by 



d t A = -E a , (62) 
7 c B d . 

where the metric operator H ab associated to the classical 3-metric q a b is 



d t {H ab E b ) = H ab e b cd V c B d 



H ab = ^fiM{q) (q ab + 2^^M(qj e abc V c ) . (63) 

In the case under consideration, the metric and the electromagnetic field 
are both given as perturbative expansions, i.e., keeping only the first non- 
trivial correction 



Qab = <5afe + e 2 7afc, (64) 

F, v = e F« + e 3 F£), (65) 

Tab 



and therefore, the metric operator H is also given in perturbative way 

H ab = H (0)ab + e 2 H (2)ab ^ 

with 



H {o)ab = 6 ab + 2 £e abc d c , (67) 

H mab = Q 5 o& 7 _ ^ + ^ ^~abc + + fP^f) . (68) 

In these expressions 7 is the trace of 7 a ft, i.e., 7 = <5 af> 7 a fc, and e abc and 

(2) h 

Va are the first non trivial corrections to e and the covariant derivative 
respectively, 
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(69) 



while Va applied to a co- vector C b is given by 



vi 2) a = ~{d albd + d blad - d dlab )5 cd c c . 



(70) 



Notice that the zeroth order of the metric operator (eq. fl67|) ) is, as ex- 
pected, just the flat operator obtained in previous works [61 [T2] . 

On the other hand, the perturbative expansion (RM|) gives rise to similar ex- 
pansions for the electric and magnetic fields (see eqs. (??) and (I55p ) which, 
inserted in the field equations (1621) lead to 



for the first order, corresponding to the free Maxwell field propagating in a 
flat background (here we have adopted vectorial notation for simplicity), and 



dtiH^E^) - H^ ab e b cd 8 c Bf = -\d t {H^ ab E { b l) ) - H^ ab t b cd d c B { d 1] ] 

+ H^ ab [e b cd d c + e^Vf]^ , 



for the correction generated by back reaction effects. The equations for the 
free field f JTTj) and f!72|) coincide, of course, with the ones obtained in [6] and 
|12j . and preserve Lorentz invariance, while the above corrections will break 
that symmetry. The easiest way to see this is by considering the wave-like 
equation to be satisfied by the potential (eq. (TH1) ). where it becomes clear 
which term breaks covariance. From this expression we can see that Lorentz 
Invariance will be broken whenever the time derivative of det(q) is different 
form zero. In our perturbative approach this quantity is given by 



whose time derivative is in general non vanishing. Hence, even for a flat 
background, if we take into account back reaction effects on the metric, 



d t (E^ + 2£V x £«) 



V x (£? (1) + 2£V x £? (1) ), 



(71) 
(72) 



Vdet(q) = 1 + - e 2 7 



(73) 
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the resulting propagation equations for the electromagnetic field will break 
Lorentz Invariance. 

Just for completeness, this wave equation in the case of interest reads, for 
the free Maxwell field 

= 0, (74) 
and, for the back reaction correction, 

[(1 -2£Vx)DA (3) ] a = -d t {H^ ab EP) 

+ {e b cd {H^ ab Vf + H^ ab d c ) + e b cd H^ ab d c ) . 

4 Summary and conclusions 

We have studied the propagation of light in two different scenarios 

1. On an arbitrary quantum curved spacetime in the semiclassical approx- 
imation of Loop Quantum Gravity. 

2. On a deviation from the quantum flat metric were the non trivial part 
is the back reaction effect of the Maxwell field. 

For the first part we obtained the effective interaction Hamiltonian for 
the gravitational and electromagnetic fields and derived the corresponding 
field equations, which can be combined to obtain a wave like equation for the 
Maxwell potential. 

In the particular case of a flat background, this wave equation reduces 
to the usual Lorentz invariant propagation. This result is also valid for any 
stationary background geometry, but it is no longer true for the general 
case, in which the wave equation contains a covariance breaking term that is 
related to the time derivative of the determinant of the 3-dimensional metric. 

As an example of this we studied light propagation in flat FRW cosmol- 
ogy dominated by matter, and solved the wave equation to obtain an effect 
that in principle can be observed, namely, that the polarization direction of 
an initial linearly polarized plane wave rotates with a frequency dependent 
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angle. However, it is not clear if the assumptions made to obtain these re- 
sults were too restrictive. It could happen that the assumptions cannot be 
maintained for the time of flight of the photons and thus there are no physi- 
cal predictions to be made. On the other hand, if the set of assumptions are 
valid there are observational consequences, such us the loss of polarization of 
a linearly polarized wave packet with a frequency spectrum. Note also that 
the polarization direction has a linear dependence on the photon energy, a 
different result from that obtained by Gambini and Pullin where the depen- 
dence is quadratic pj. However, the fact that we do observe light with a 
large amount of linear polarization tells us that, if this effect actually exists, 
it must be much smaller than expected, and, moreover, using recent observa- 
tional data it is possible to put a very severe bound on the phenomenological 
constant £ [3]. 

A second and for us more important problem was to analyze the propa- 
gation of light taking into account back reaction effects. We have seen that 
Lorentz invariance is also broken and, although we have not solved the equa- 
tions, it is reasonable to believe that wave propagation would present similar 
effects to the ones obtained for a FRW background. However, in this case, 
any induced effect, if present, would be much more difficult to observe since 
it is of higher order: the corrections are second order in the small parameter 
e and, besides, of order ££ p . 

It is worth mentioning that polarized light travelling on a media with 
an index of refraction induced by the quantum spacetime is a very sensitive 
tool to study these corrections and it is a worthwhile problem to obtain a 
predicted value for the rotation of the polarization direction. 

As a final comment we would like to mention that it is surprising to ob- 
tain noncovariant semiclassical equations of motion coming from a covariant 
formalism. Loop quantum gravity is by construction a covariant theory, al- 
though the 3+1 splitting hides this fact. The only possible place where this 
covariance can be broken is in the use of semiclassical states. These states 
are only gauge invariant with respect to the rotation group on the spacial 
surface but do not satisfy the hamiltonian constraint (otherwise they cannot 
be peaked around a classical metric). Maybe it is impossible to define a 
covariant semiclassical approximation of loop quantum but this would be a 
rather unwanted feature on the theory. 
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